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1.  INTRODUCTION  -•* 


such  diverse  phenomena  as 


turbulence,  extreme  properties  of  stochastic  systems,  the  statistics  of  nonlinear 
wave-wave  interactions  and  the  detection  of  clustered  events.  One  of  the 
theories  developed  to  explain  the  transition  from  laminar  to  turbulent  flow  is 
through  the  use  of  intermittent  turbulent  bursts  interacting  with  the  laminar 
flow  to  generate  more  bursts  and  thus  form  a  cascade  to  a  fully  developed  tur¬ 
bulent  state.  Thus,  the  statistical  distribution  of  such  bursts  is  important  in 
determining  the  rate  of  transition  from  laminar  to  turbulent  fluid  flow  and  can 
thereby  strongly  influence  the  drag  properties  of  fluid  flowing  past  a  body.  A 
second  process  in  which  intermittency  is  important  is  the  anomolous  transport 
of  charge  in  amorphous  materials.  The  current  flow  in  these  materials  exhibits 
a  clustering  of  events  in  time  not  unlike  that  observed  in  fluid  turbulence.  In 
both  of  these  systems  the  physical  observables,  e  g.  the  correlation  function 
exhibit  scaling  behavior  with  exponents  characteristic  of  the  process.  A  Anal 
example  is  the  stress  relaxation  of  polymers,  in  which  case  the  scaling  behavior 
effects  such  material  properties  as  brittleness,  dielectric  loss,  etc.  -Under  this  ^ 
.-contract  we  have  developed  a  new  random  walk  model^rhich  greatly  increases 
our  understanding  of  such  processes. 


Many  scaling  relations  for  complex  systems  in  the  physical  sciences  involve 
non-integer  exponents.  We  list  several  examples  of  these  below,  which  although 
well  known  have  provided  us  with  the  orientation  necessary  to  investigate  more 
unfamiliar  physical  situations.  We  interpret  the  notHnteger  exponents  in  the 
scaling  relations  as  indications  of  singularities  arising  from  probability  distribu¬ 
tions  with  long  tails  governing  the  physical  observables.  The  existence  of  such 
distributions  to  describe  turbulence,  electron  transport  in  amorphous  materials 
and  other  areas  has  been  proposed  by  a  number  of  investigators.  But 


-1- 


-2- 


heretofore  no  systematic  investigation  of  the  physical  implications  of  these  dis¬ 
tributions  has  been  presented.  The  initial  work  in  this  area  has  been  done  under 
this  contract  in  the  sequence  of  fifteen  papers  described  below. 

If  the  first  moment  of  the  appropriate  physical  observable  with  respect  to 
one  of  these  long  tailed  distributions  diverges,  then  no  scale  exists  in  which  to 
gauge  measurements  and  structure  occurs  on  all  scales.  The  concepts  of  self¬ 
similar  fractals,  non-dillerentiability  ,  and  also  noninteger  exponents  all  accom¬ 
pany  the  divergence  of  such  low  order  moments.  Ihe  analysii  we  have  done 
centers  around  the  construction  of  simple  mechanical  modeli  exhibiting  these 
properties.  In  particular  we  examine  random  walkexamples  where  the  above 
characteristics  appear  simply  and  naturally.  Thesi  random  processes  have  an 
inherent  self-similar  (fractal)  scaling  in  space,  timi,  frequency  or  other 
appropriate  variable.  They  can  be  used  to  model  complex  systems  of  interest 
which  exhibit  features  spanning  many  decades  of  stale.  As  mathematicians  we 
categorize  possible  behaviors  and  as  physicists  we  tearch  for  the  basic  interac¬ 
tions  which  lead  to  such  behavior. 
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2.  SURPRISES  IN  EXPERIMENTAL  PHYSICS 

When  one  begins  the  study  of  physics,  one  invariably  discovers  simple  laws 
such  as  momentum  (mv  =  constant)  conservation  and  energy  conservation 

(  h-mv  2  =  constant )  .  One  also  encounters  many  well-defined  concepts,  such 

as  the  mobility  /*,($  =  jl  •  2),  say  for  a  conductance  electron  in  a  metal.  Sim¬ 
ple  relations  are  ubiquitous,  e.g.,  the  mean  square  displacement  of  a  random 
walker  after  N  steps  is  proportional  to  N;  the  vibrational  density  of  states  p(u) 
of  a  D  dimensional  crystal  at  low  frequencies  scales  as  uD~l,  etc.  Everywhere 
exponents  seem  to  be  integers  and  quantities  are  well-defined. 

When  one  deals  with  actual  systems  which  are  sufficiently  complex,  it  is 
often  the  case  that  the  best  description  is  in  terms  of  probability  distributions. 

A  Poisson  distribution  is  characterized  by  its  first  moment,  and  a  Gaussian  dis¬ 
tribution  by  its  first  two  moments.  When  distribution  have  long  enough  tails, 
the  first  few  moments  will  not  characterize  the  distribution.  Distributions  with 
infinite  moments  will  lead  to  physics  described  by  non-integer  exponents,  and  to 
surprises  which  run  counter  to  our  intuition.  We  list  some  examples  which  we 
have  investigated  where  non-integer  exponents  abound. 

A.  High  Velocity  Impacts1'^ 

An  important  materials  science  problem  is  the  effect  of  high  velocity  micro 
meteorites  (mass  m,  velocity  v)  on  spacecraft.  As  a  result  of  a  collision  many 
interesting  phenomena  occur,  deformation  of  the  meteor,  (projectile)  crater 
formation  on  the  surface  of  the  spacecraft  (target),  shock  wave  generation, 
melting,  and  crack  propagation.  One  would  suspect  that  the  crater  formed  by 

the  collision  would  be  proportional  to  the  kinetic  energy  of  the  meteor  ( —mu2) 

leading  to  the  radius  of  the  crater  R  scaling  as  R8  ~  v  2  Momentum  conserva¬ 
tion  is  not  considered  because  target  material  is  thrown  backwards  during  the 
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crater  formation  thus  enhancing  the  forward  momentum.  The  surprise  is  that 
the  crater  depth  scales  as  a  power  of  the  velocity,  i.e. 

R  ~va  .  a  =  0.56  (1) 

but  a  *  2/3.  This  implies  a  new  conservation  law  given  by 

mu1-7*  =  constant  (2) 

Ve  have  been  able  to  explain  this  scaling  behavior  by  using  a  random  walk  model 

for  the  movement  of  dislocations  in  the  target  material. 

&  Charge  Transport  in  Amorphous  F  ‘ns3-5 

In  the  xerographic  process  a  flash  of  light  creates  a  layer  of  electron-hole 
pairs  in  an  amorphous  film.  The  charges  separate  under  the  influence  of  an 
external  electronic  Held  E.  If  the  film  is  of  depth  L,  and  there  exists  a  mean 
transit  time  T  for  charges  to  cross  the  film,  then  one  expects  that  L/T  =  juE. 
Instead  one  finds  (for  a-As2  Seg  )  that 

T~l  ~  (E/L)22  =  ( E/L)l/A0  (3) 

for  a  wide  variety  of  E  and  L  values.  The  current  l(t)  generated  by  the 

charge  movement  is  found  to  have  two  regimes.  In  the  early  phase 

/(t)~f“‘*«  (4a) 

and  in  the  late  phase  (caused  by  absorption  at  the  far  surface) 

/(f)~f"l"4B  (4b) 

A  continuous  time  random  walk  model  is  the  only  complete  description  of  this 

1  O 

phenomenon. 

G.  Polymers6 

1.  The  world  of  polymers  is  hill  of  fractional  exponents.  The  mean  square  dis¬ 
placement  <  R2  (N)  >  of  a  self-avoiding  random  walk  (SAW)  after  N  steps 
can  be  used  to  model  the  end-to-end  distance  of  a  polymer,  in  a  good  sol¬ 
vent.  In  three  dimensions  (3D) 
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<  R*(N)  >~fP  (5) 

y 

with  y  ~  1.2.  Mandelbrot  has  suggested  that  y  =  2/p  where  p  is  the 

Q 

’'fractal'’  dimension  of  the  SAW  and  is  approximately  In  6/ln3  ~  1.6. 

o  21 

Shlesinger  has  used  some  recent  results  of  Seshadri  and  West  to  explain 
this  phenomenon  using  the  newly  developed  concept  of  fractal  random 
walks  developed  under  this  contract®’  12,15.18 

2.  The  vibrational  density  of  states  p(o)  of  certain  proteins  is  found  at  low  fre¬ 
quencies  to  be® 

p(o)~uMt0*  (0) 

This  is  consistent  with  the  protein  having  a  3D  SAW  shape  and  p(u)  ~  uP'  *. 

and  in  contrast  to  the  Debye  law  p{u)  ~  uP~l  =  u2  in  3D. 

3.  The  repetition  time  tr  for  polymer  melts  is  defined  as  the  cross-over  time 
from  rubber-like  to  liquid-like  behavior.  It  is  found  to  scale  with  polymer 
mass  M  as 

rR  "*  HI™  (7) 

O 

A  simple  diffusion  model  for  the  motion  of  constrained  polymer  chains 
would  yield  Tr  ~M2. 

4.  The  dielectric  response  of  polymers  (as  well  as  many  other  materials)  is 
governed  by  non-integer  exponents.  The  frequency  dependent  dielectric 
constant  e  (  cj)  is  given  by 

(8> 

£m  *  t0  .N  ® 

Williams  and  Watt10  find  that  a  good  form  for  the  response  function  is 

?(f)  =  M  (XosBl  (9) 

The  relation  between  this  response  function  and  long  tail  distributions  has 

1 A 

been  discussed  In  some  detail  by  Montroll  and  Shlesinger. 
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5.  Amorphous  solids  are  not  in  thermodynamic  equilibrium  at  temperatures 
below  their  glass  transition  temperature.  These  non-equilibrium  states  are 
metastable  and  slowly  relax  towards  equilibrium  and  in  so  doing  affect 
many  properties  of  the  material,  e.g.  the  material  becomes  stiller  and  more 
brittle,  its  damping  decreases  and  so  do  the  creep-  and  stress-relaxation 
rates,  the  dielectric  constant,  the  dielectric  loss  etc.  The  lifetime  of  such 
metastable  states  when  the  transition  is  induced  by  correlated  fluctuations 
has  been  studied  by  West  and  Linde nberg. 19,20  The  probability  density  in 
this  case  is  shown  to  satisfy  a  Fokker- Planck  equation  having  a  diffusion 
coefficient  which  is  both  time  dependent  and  dependent  on  the  state  of  the 
system. 
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a  FROM  NONtWTKGER  EXPONENTS  TO  SELF-SIMILARnT 

Integer  exponents  can  usually  be  traced  back  to  the  analytic  behavior  of  an 
appropriate  function  which  can  be  expanded  in  a  Taylor  series.  Non-integer 
exponents  imply  the  presence  of  singularities  and  the  breakdown  of  a  Taylor 
series  due  to  the  divergence  of  a  coefficient.  Our  main  theme  is  that  singulari¬ 
ties  and  thus  non-integer  exponents  arise  in  complex  systems  because  they  exhi¬ 
bit  randomness  on  many  scales.  As  an  example,  consider  the  probability 
(t)  dt  for  an  event  (electron  hopping,  dislocation  movement,  crack  propaga¬ 
tion,  etc.)  to  occur  in  the  time  interval  (t,  t+dt).  Its  Laplace  transform  is  defined 
by 

m 

it*(s)  *  f  «“**  i/(t )dt  (10) 

o 

If  all  it*  moments  <fB>  of  ii  (t)  exist,  then 

*  i  +  £  *-£-<*">  (id 

n*l  n! 

and  times  will  be  measured  in  units  of  <t>.  If  however.  <t>  is  infinite,  i.e., 

<t>  =  £  t  i>(t)dt  =  -  *  (s  =0)  =  - 

then  t/<t>  is  not  an  appropriate  dimensionless  time,  and  i>*(s)  cannot  be 

represented  in  a  Taylor  series  about  s  =  0.  This  will  occur  if  at  long  times 

^r(t)  ~  t'1'*  ,  0  <  o  <  1.  A  proper  small  s  expansion2*5,11  of  ^*(s)  will  include 

an  s°  term.  Effectively  t  is  replaced  by  f "  in  all  dimensional  analysis.  This  is 

the  underlying  cause  for  the  non-integer  exponents  in  (3)  and  (4).  Such  (t) 

arise  from  a  distribution  of  deep  traps  in  the  material  which  can  capture 

o  93 

charges  in  amorphous  materials,  dislocations  in  elastic  materials  '  and/or 
excitations  in  polymers.®8,27  There  Is  no  average  time  in  these  problems  when 
<t>  =  »,  so  that  the  duration  between  events  occurs  on  all  time  scales.  A  self¬ 
similar  set  of  burst  of  events  has  this  property. 


Systems  with  inherent  scaling  are  most  naturally  described  by  scaling 
equations,  and  not  by  differential  equations.  In  fact,  the  appropriate  functions 
may  be  nowhere  differentiable.  Scaling  brings  us  into  the  province  of  the  renor¬ 
malization  group  (RNG). 1 1  Let  us  consider  how  the  free  energy  F  of  a  system 
scales  near  a  critical  point. 

Let  K  describe  interactions  within  a  small  system,  and  K'(K)  interactions 
between  these  small  systems.  When  a  non-trivial  fixed  point  of  the  RNG  transfor¬ 
mation  1C(K)  =  K  exists,  one  may  switch  to  a  scaling  variable,  u.  to  find 

F(u)  =  rDF(\u)  +  G(u)  (12) 

where  in  D  dimensions,  u  describes  the  small  system  of  unit  size,  X  u  the 

larger  system  of  size  lD.  X  is  a  relevant  eigenvalue  of  the  RNG  transformation, 

and  G  is  an  analytic  function.  Let  us  iterate  (12)  n  times  and  then 

let  n  -»  «.  to  obtain 

F{u)  =  lim  t^  F(\nu)  +  y?  l-i°  G{\*u)  .  (13) 

l  i-o 

It  is  usually  assumed  that  the  fist  term  on  the  rhs  vanishes  as  n  ,  and  thus 
all  singularities  in  F(u)  must  shift  to  the  term  involving  the  summation  over  G.  A 
solution  to  the  singular  part  of  (13)  is 

F{u)  =  A(u)lu  (14a) 

where  A(u)  is  oscillatory  in  In  u  with  period  In  A,  i.e. 

A(u)  =  A(Au)  =  £*,  (14b) 

« 

In  our  work  we  introduced  fractal  random  walks1 1-15  processes  which  produce 
non-integer  exponents  naturally  in  the  form  of  fractal  dimensions,  and  integral 
transforms  of  probability  distributions  exhibit  precisely  the  scaling  of  (12)  - 
(14). 

We  treat  separately  the  most  famous  scale-invariant  problem  the  "l/f 
noise”  problem.1®'17' 10  For  many  diver  e  systems  the  power  spectrum  S(f) 


(the  Fourier  transform  of  a  second  order  correlation  function  C  (t) )  has  a  low 
frequency  component  scaling  as  the  inverse  of  the  frequency,  i.e. 

ee 

S(f)  =  Ref  e9*^1  c(t)dt  ~  1/ /  os  /  ->0 
0 

In  this  regime  the  integrated  power  spectrum  is  independent  of  scale.  i.e. 

fT  il  -  fT  <*(///.) 

/L  f  "  /L  r/f. 

For  a  completely  random  system  C(r)  =  a~l/r  where  r  is  a  relaxation  time.  In 
a  complex  system  a  distribution  of  relaxation  times  p(r)  can  exist.  Then 

5(/)  =  f  e~t/r p(r)dr  . 

0 

If  p(r)  ~  t~1  then  a  regime  of  1/f  noise  will  result.  We  show  that  such  a  p(t) 
arises  when  the  probability  of  fluctuation  depends  on  the  product  of  several  ran¬ 
dom  variables.16  In  this  case  p(r)  will  be  a  log-normal  distribution  in  contrast 
to  the  normal  distribution  which  arises  for  the  distribution  of  a  sum  of  random 
variables.  Our  p(r)  involves  a  geometric  mean  while  the  normal  distribution 
involves  a  arithmetic  mean. 

The  log  normal  distribution  also  describes  (at  longer  times)  the  lifetime  of 
many  materials.  We  are  at  present  investigating  deviations  to  this  law  at  early 
times. 

The  random  walk  models  developed  to  understand  the  phenomenon  of  clus¬ 
tering  in  either  space  and/or  time  have  relied  on  long  range  transition  proba¬ 
bilities  for  stepping  between  sites  on  a  lattice.  In  the  continuum  limit  the  proba¬ 
bility  density  has  been  shown  to  satisfy  an  integral-differential  equation  rather 
than  a  Fokker-Planck  equation.  *2. 14, 15  Tjae  80iuti0-j  t0  this  equation  is  a  Ldvy 
distribution  having  the  characteristic  function  9»(k)  =  e  ")’iklMi.e.  the  Fourier 
transform  of  the  probability  density.  The  Ldvy  distribution  describes  intermit¬ 
tent  or  clustered  processes  having  mean  first  passage  times  and  maxima 
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moments  which  scale  with the  characteristic  exponent  \i.  The  Ldvy 

parameter  p.  has  been  identified  with  the  fractal  dimensionality  for  the  pro- 

cess  1  and  can  be  expressed  in  terms  of  the  random  walk  parameters.  ’  This 

2  23 

Identification  has  been  us«d  in  both  the  high-velocity  impact  problem  '  and 

B 

the  end-to-end  distance  of  a  polymer  chain  . 

The  dynamic  equatioas  describing  the  evolution  of  such  systems  are  neces¬ 
sarily  nonlinear  and  stochastic,  as  for  example  the  description  of  the  generation 
and  evolution  of  dislocations  in  elastic  materials.^4  The  analysis  of  such  non¬ 
linear  stochastic  equations  is  notoriously  difficult  so  we  have  also  developed 
under  this  contract  a  linesrization  method  which  is  instantaneously  optimal  at 
all  stages  of  evolution  of  tie  process.  Detailed  comparison  of  the  results  of  this 
technique  with  Monte  Carle  calculations  on  certain  model  systems  are  presented 

pc 

by  West.  et.  al.  and  the  comparison  is  quite  favorable. 

Energy  transport  in  nolecular  aggregates,  e.g.  polymeric  materials,  has 
long  been  a  useful  diagnostic  for  the  determination  of  their  structural  and 
dynamical  properties.  On*  of  the  serious  limitations  of  the  previous  theories  has 
been  the  lack  of  understaiding  of  the  temperature  dependence  of  the  measured 
transport  properties  of  th*  aggregates.  The  limitation  arises  precisely  because 
these  have  been  infinite  temperature  theories.  West  and  Lindenberg^®’^7  have 
developed  a  novel  approach  for  the  incorporation  of  finite  temperature  effects 
into  a  phenomenologically  based  model.  By  including  the  appropriate  dissipa¬ 
tive  effects  in  the  equations  for  exciton  dynamics  in  molecular  aggregates,  the 
new  theory  is  correct  at  finite  (arbitrarily  low)  temperatures  and  indications  are 
that  it  will  be  able  to  expUin  a  number  of  "anomolous"  temperature  dependen¬ 
cies  of  spectral  and  transfort  properties. 
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